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Abstract. A systematic low-energy effective field theory for hole-doped antiferromagnets on 
the honeycomb lattice is constructed. The formalism is then used to investigate spiral phases 
in the staggered magnetization as well as the formation of two-hole bound states. 
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1. Motivation 

Although the phenomenon of high-temperature superconductivity was discovered more than 
twenty- five years ago [1] , the dynamical mechanism behind it still remains a mystery today. The 
reason may be readily identified: Due to the nonperturbative nature of the problem, analytic 
studies usually suffer from uncontrolled approximations, while numerical simulations of the 
microscopic Hubbard and t- J- type models suffer from a severe fermion sign problem away from 
half-filling. 

In the present work, we will address the physics of the antiferromagnetic precursors of high- 
temperature superconductivity from a universal and model-independent point of view, based on 
the method of effective Lagrangians. The effective Lagrangian method is very well-established in 
particle physics, where chiral perturbation theory (CHPT) [2j [3] represents the effective theory 
of quantum chromodynamics (QCD). It operates in the Goldstone boson sector of QCD, taking 
into account the octet of the pseudoscalar mesons that emerge due to the spontaneously broken 
chiral symmetry. Baryon chiral perturbation theory |H [5j EJ [7] is an extension of CHPT as it 
also includes heavy particles such as the baryon octet or decuplet, i.e. those degrees of freedom 
that remain massive in the chiral limit. 

Spontaneous symmetry breaking is also a very common phenomenon in condensed matter 
physics: Magnons e.g. are the Goldstone bosons of a spontaneously broken spin symmetry 
SU(2) S — > U(l) s , while phonons emerge due to a spontaneously broken translation symmetry. 



For these systems effective Lagrangians have been constructed in Refs. [8J [9l [TU1 [TTj and various 
applications have been considered in Refs. [12j [131 EH EH IT7] ITS]. 

However, the incorporation of heavy degrees of freedom in the condensed matter domain 
- in a transparent and systematic manner - has only been achieved recently. An effective 
theory for weakly doped antiferromagnets, the analog of baryon chiral perturbation theory, 
has been constructed for hole- and electron-doped antiferromagnets on the square lattice in 
Refs.plJJ l20l |2T] and for hole-doped antiferromagnets on the honeycomb lattice in Ref. |22j. 
Further applications of the formalism were discussed in Refs. [231 12^1 125] . demonstrating that 
the effective Lagrangian method allows one to gain insight into the physics of these insulating 
precursors of high-temperature superconductors in a systematic and unambiguous manner. 

In this presentation we will focus on hole-doped antiferromagnets on a honeycomb lattice. 
We first review the systematic construction of the effective Lagrangian for magnons and holes 
for these systems and will then consider two applications of the formalism: The emergence of 
spiral phases in the staggered magnetization order parameter as well as the formation of two-hole 
bound states mediated by magnon exchange. 

We would also like to point out that in a recent article on an analytically solvable 
microscopic model for a hole-doped ferromagnet in 1+1 dimensions [26], the correctness of the 
effective field theory approach was demonstrated by comparing the effective theory predictions 
with the microscopic calculation. Likewise, in a series of high-accuracy investigations of 
the antiferromagnetic spin-| quantum Heisenberg model on a square lattice using the loop- 
cluster algorithm |27| |2"E1 1251 130] . the Monte Carlo data were confronted with the analytic 
predictions of magnon chiral perturbation theory and the low-energy constants were extracted 
with permille accuracy. All these tests unambiguously demonstrate that the effective Lagrangian 
approach provides a rigorous and systematic derivative expansion for both ferromagnetic and 
antiferromagnetic systems. 

2. Construction of the effective field theory for holes and magnons 

The effective Lagrangian approach is based on a symmetry analysis of the underlying theory, 
i.e., in the present case the Hubbard or t-J-type models which are believed to be minimal models 
for high-temperature superconductors [31] . We thus first perform a rigorous symmetry analysis 
of these models in order to then construct the effective Lagrangian for holes and magnons. 

2.1. Symmetry analysis 

Let c xs denote the operator which creates a fermion with spin s E {7, 4} on a lattice site 
x = (a?x,X2). The corresponding annihilation operator is c xs . These fermion operators obey the 
canonical anticommutation relations 

{ c xsi c ys'} = $xy$ss'i { c xs: Cys'} = { c xsi C ys'} = ^- (1) 

The second quantized Hubbard Hamiltonian is defined by 

H = —t ^2 { c xs c ys + Cys c xs) + C x J [ Cx t C xi C xi ~ I 1 ^2 C xs c xs, (2) 

where {x, y) indicates summation over nearest neighbors, t is the hopping parameter, and the 
parameter U > fixes the strength of the Coulomb repulsion between two fermions located on 
the same lattice site. The parameter y! denotes the chemical potential. 

Using the fermion creation and annihilation operators, we define the following SU(2) S Pauli 
spinors 

4=(4t,cy, *-(£)■ W 




Figure 1. Bipartite non-Bravais honeycomb lattice consisting of two triangular Bravais 
sublattices. The translation vectors are a± and a% . 



In terms of these operators, the Hubbard model can be reformulated as 

h = -tj2( c Uy + c\c x ) + — 53(4pas - 1) 2 - //XX4 c x ~ 

(xy) x x 



(4) 



The parameter \x = // — ^ controls doping where the fermions are counted with respect to 
half-filling. 

Since all terms in the effective Lagrangian must be invariant under all symmetries of the 
Hubbard model, a careful symmetry analysis of Eq.@ is needed. Let us divide the symmetries of 
the Hubbard model into two categories: Continuous symmetries, which are internal symmetries 
of Eq.@, and discrete symmetries, which are symmetry transformations of the underlying 
honeycomb lattice depicted in figure 1. 

The continuous symmetries are the SU(2) S spin rotation symmetry, the U(1)q fermion 
number symmetry and its non-Abelian extension SU(2)q. The discrete symmetries include 
translations (Di) along the two primitive lattice vectors of the honeycomb lattice, rotations (O) 
by 60 degrees around the center of a hexagon, and reflections (R) at the xi-axis going through 
the center of the hexagon. There is also time reversal which is implemented by an anti-unitary 
operator T. This symmetry will be considered later on in the effective field theory framework. 
Furthermore, in the construction of the effective field theory for magnons and holes, it will turn 
out to be useful to incorporate the combined symmetry O 1 consisting of a spatial rotation O 
and a global SU(2) S spin rotation g = ia%. Likewise, we define the composed transformation 
T", consisting of a regular time-reversal T and the specific spin rotation g = io<i- 

In |32} 133] . Yang and Zhang proved the existence of a non-Abelian extension of the U(1)q 
fermion number symmetry in the half-filled Hubbard model. This pseudospin symmetry contains 
U(1)q as a subgroup. The SU(2)q symmetry is realized on the square as well as on the 
honeycomb lattice and is generated by the three operators 
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Q = 2( -1 )*c*je»t> 



Q 3 = E^4 t c -t + 4^-i) = ^Q- (5) 

X 

The factor (— l) x distinguishes between the two triangular sublattices A and B of the honeycomb 
lattice. Defining Q 1 and Q 2 through = Q 1 ± iQ 2 , one readily shows that the SU(2)q Lie- 
algebra [Q a ,Q b ] = i£ a bcQ c , with a,b,c G {1,2,3}, indeed is satisfied and that [H, Q] = with 
Q = (Q , Q 2 ,Q S ) for the Hubbard Hamiltonian with /i = 0. 

In order to write the Hubbard Hamiltonian Eq.([2|) or Eq.@ in a manifestly invariant form 
under SU(2) S x SU(2)q, we arrange the fermion operators in a 2 x 2 matrix-valued operator, 
arriving at the fermion representation 

V «i4 -(- 1 ) 4 t / 

This allows us to write down the Hubbard Hamiltonian in the manifestly SU(2) S , U(1)q, Di, 
O, O' and R invariant form 

H = -\Y, nC x C x+i + C[ + f x ] + ^ £ nctC x ClC x ] -f£ TrY :,/>,;. (7) 

x,i x x 

The (T3 Pauli matrix in the chemical potential term prevents the Hubbard Hamiltonian from 
being invariant under SU(2)q away from half-filling. For \i ^ 0, SU(2)q is explicitly broken 
down to its subgroup U(1)q. In addition, the pseudospin symmetry is realized in Eq.([7]) only for 
nearest-neighbor hopping. As soon as next-to-nearest-neighbor hopping is included, the SU{2)q 
invariance gets lost even for fj, = 0. The continuous SU(2)q symmetry contains a discrete 
particle-hole symmetry. Although this pseudospin symmetry is not present in real materials, it 
will play an important role in the construction of the effective field theory. The identification 
of the final effective fields for holes will lead us to explicitly break the SU(2)q symmetry in 
subsection \~ 



2.2. Effective field theory for magnons 

In this section we review the construction of the effective theory for antiferromagnetic magnons 
on the honeycomb lattice. The basic object in the effective theory is the staggered magnetization 
order parameter of the antiferromagnet, which is described by a unit-vector field 

e(x) = (sm0(x) cos <p(x), sin 9{x) sin</?(x), cos 9(x)), (8) 

in the coset space SU(2) S /U (l) s = S 2 , where x = (x\,X2,f) denotes a point in (2+l)-dimensional 
space-time. A key ingredient for constructing the effective field theory is the nonlinear realization 
of the global SU(2) S spin symmetry which is spontaneously broken down to its U(l) s subgroup, 
which is discussed in detail in Ref. [19]. 

It turns out that the above parametrization of the magnon degrees of freedom through the 
vector e(x) is not appropriate for this construction. Rather, one uses composite magnon fields 
Vfj,(x) whose components will be used to couple the magnons to the fermions. The composite 
magnon field is defined by 

v^x) = u(x)d fl u(x) J< , (9) 
where the matrix u(x) is related to the original magnetization vector e{x) = (e\(x), e2(x),es(x)) 
by 

l + e 3 (x) ei{x) - ie 2 (x) 



u(x) 



y/2(i+e 3 (x)) \ -ei(x) - ie 2 (x) l + e 3 (x) 

cos sin (^f^J exp(—iip(x)) 

(*M) «pM*)) cos (¥) 



sin 



(10) 



The coupling of magnons to holes is then realized through the matrix-valued anti-Hermitean 
field 

Vpix) = iv^(x)a a , v±(x) = v\{x) =F iiS^x), (11) 

which decomposes into an Abelian "gauge" field v^(x) and two vector fields v^(x) "charged" 
under the unbroken subgroup U(l) s . Here a are the Pauli matrices. These fields have a well- 
defined transformation behavior under the symmetries which the effective theory inherits from 
the underlying microscopic t-J model: 



SU(2) S 
Di 
O 



R 
T 



Vp(x)' = h(x)(v^(x) + <9 M )/i(x) t , 
Di v^(x) = Vfj,(x), 



it 



° Vl (x) = t{Ox)[\ Vi {Ox) + &v 2 {Ox) + \d x + ^<9 2 }t(Ox 
°v 2 (x) = t(Ox){ - ^v x (Ox) + \v 2 {Ox) - &d x + i5 2 }r(Ox) t , 

°v t {x) = T(Ox)(v t (Ox) + d t )T{Ox)\ 
R vi(x) = Vi(Rx), R v 2 (x) = -v 2 (Rx), 
R v t (x) = v t {Rx), 

T Vi(x) = T(Tx)(vi(Tx) + di)T{Tx)\ 



T v t {x) = -T{Tx){v t {Tx) + d t )r(T x y. 
In the above symmetry transformations, we have introduced the matrix r(x), 



t(x) 



— exp(— i<p(x)) 

exp(i(p(x)) 



(12) 



(13) 



Finally, the Abelian "gauge" transformation 

h{x) = exp(ia(x)o"3) (14) 
belongs to the unbroken U(l) s subgroup of SU(2) S and acts on the composite vector fields as 

"2 s )' = uj(x) -dfta(x), 



v^{x) exp(±2ia(x)). 



(15) 



The magnon action, in terms of the composite magnon field v^(x), can now be expressed as 



S[v* 



J d 2 x dt 2p s ^vfv i + X 



2 V i v t 



(16) 



Although the expression v^v looks like a mass term of a charged vector field, it is just the 
kinetic term of a massless Goldstone boson, since it contains derivatives acting on u{x). 



2.3. Effective field theory for holes and magnons 

Analytic calculations as well as Monte Carlo simulations in i-J-like models on the honeycomb 
lattice have revealed that at small doping holes occur in pockets centered at lattice momenta 

f,a = _ k f3 _ ( Q) and their 

copies in the periodic Brillouin zone \34\ 135]. The honeycomb 
lattice, illustrated in figure 1, is a bipartite non-Bravais lattice which consists of two triangular 
Bravais sublattices. The corresponding Brillouin zone and the location of the corresponding 
hole pockets are shown in figure 2. The single-hole dispersion relation for the t-J model on the 
honeycomb lattice is illustrated in figure 3. 




Figure 3. Energy-momentum dispersion relation Eh(k)/t for a single hole in the t-J model on 
the honeycomb lattice for J/t = 2. 



The effective field theory is defined in the space-time continuum and the holes are described 
by Grassmann-valued fields ipl( x ) carrying a "flavor" index / = a, (3 that characterizes the 
corresponding hole pocket. The index s = ± denotes spin parallel (+) or antiparallel (— ) to the 
local staggered magnetization. Under the various symmetry operations the hole fields transform 

as 



517(2), 
U(1)q 
A 
O 



R : 



ip±(x)' = exp(±ia(x))ip±(x) , 
Di ^ ± {x) = exp(ik{ ai)ip±{x), 



°^(x) = tzMtMOx) ±i^)^(Ox), 
( T ) _ m rxdC mun( Ox) t 



=F exp(=fiip(Ox) =F i^-)i^{Ox) 



R 



^Ux) = ^i(Rx), M(^) = mRx), 



T : 



T V>±(x) = exp(Ti<p(Tx))i${Tx), 
T ip±(x) = exp(^fi(p(Tx))i/j^ (Tx), 



x) 
x) 



exp(±i(^(Tx))V ; ± (Tx), 
exp(±i(p(Tx))i/>%.(Tx). 



(17) 



Here J7(1)q is the fermion number symmetry of the holes. Interestingly, in the effective 
continuum theory the location of holes in lattice momentum space manifests itself as a " charge" 
k\ under the displacement symmetry D{. 

Now that the relevant low-energy degrees of freedom have been identified and the 
transformation rules of the corresponding fields have been worked out, the construction of the 
effective action is uniquely determined. The low-energy effective action of magnons and holes is 
constructed as a derivative expansion. At low energies, terms with a small number of derivatives 
dominate the dynamics. Since the holes are heavy nonrelativistic fermions, one time-derivative 
counts like two spatial derivatives. Here we limit ourselves to terms with at most one temporal 
or two spatial derivatives. One then constructs all terms consistent with the symmetries listed 
above. The effective action can be written as 



d 2 x dt ^2C 1lit , 



(18) 



where denotes the number of fermion fields that the various terms contain. The leading 
terms in the pure magnon sector take the form 



(19) 



The leading terms with two fermion fields (containing at most one temporal or two spatial 
derivatives), describing the propagation of holes as well as their couplings to magnons, are given 
by 



E 

f=a,fi\ s=+,- 



+iK[(D 1 + isa f D 2 )ipK{v( + isa f v s 2 )^L s 
-(uf + isa f v s 2 )i>^{D 1 + isa f D 2 )ip f _ ; 
+<TfL1>pe ij f&f + N^fjvfv-^f 

+isa f N 2 {^vlv 2 s ^ s - il> f Jv&i a il){) 



f 



(20) 



Here M is the rest mass and M' is the kinetic mass of a hole, A and K are hole-one-magnon 
couplings, while L, N±, and N 2 & re hole-two-magnon couplings. Note that all low-energy 
constants are real-valued. The sign is + for a and — for (3. We have introduced the field 
strength tensor of the composite Abelian "gauge" field 



f? j (x) = d i v*(x)-d j vf(x), 
and the covariant derivatives D t and Dj acting on ip± (x) as 



(21) 



AV4 {x 



d t ± ivt(x) - \i ipi(x), 



D^{(x) = \di±ivf(x)]fpi(x). 



(22) 



The chemical potential [i enters the covariant time-derivative like an imaginary constant vector 
potential for the fermion number symmetry U(1)q. It is remarkable that the term proportional to 
A with just a single (uncontracted) spatial derivative satisfies all symmetries. Due to the small 
number of derivatives it contains, this term dominates the low-energy dynamics of a lightly 
hole-doped antiferromagnet on the honeycomb lattice. Interestingly, for antiferromagnets on 
the square lattice, a corresponding term, which was first identified by Shraiman and Siggia, 
is also present in the hole-doped case |20j . On the other hand, a similar term is forbidden 
by symmetry reasons in the electron-doped case [21 J. For the honeycomb geometry we even 
identify a second hole-one-magnon coupling, K, whose contribution, however, is sub-leading. 
Interestingly, hole- or electron-doped antiferromagnets on the square lattice do not allow terms 
containing the field-strength tensor /y in £2- 

Finally, the leading terms without derivatives and with four fermion fields are given by 

*=+,- 

+ G 2 CV>f $ + G 3 CV^V^}- (23) 

The low-energy four- fermion coupling constants G%, G2, and G3 are also real- valued. Although 
potentially invariant under all symmetries, terms with two identical hole fields vanish due to the 
Pauli principle. 

It is important to note that in the above construction of the effective Lagrangian for the 
t-J model which contains holes only, a crucial step was to identify the degrees of freedom that 
correspond to the holes. In order to remove the electron degrees of freedom one has to explicitly 
break the particle-hole SU(2)q symmetry, leaving the ordinary fermion number symmetry U (1)q 
intact. This task can be achieved by constructing all possible fermionic mass terms that are 
invariant under the various symmetries. Picking the eigenvectors which correspond to the lowest 
eigenvalues of the mass matrices then allows one to separate electrons from holes. This procedure 
is described in detail in Ref. [22J 

We would like to point out that the leading order terms in the effective Lagrangian for 
magnons and holes constructed above exhibit two accidental global symmetries. First, we notice 
that for c — > 00 and without the term proportional to iK in £2, Eqs. (ll9l) and (I20p . have an 
accidental Galilean boost symmetry. Although the Galilean boost symmetry is explicitly broken 
at higher orders of the derivative expansion, this symmetry has physical implications: The 
leading one-magnon exchange between two holes, to be discussed in the next section, can be 
investigated in their rest frame without loss of generality. 

In addition, we notice an accidental global rotation symmetry 0(7). Except for the term 
proportional to iK, £2 of Eq. (l20j) is invariant under a continuous spatial rotation by an angle 
7. This symmetry is not present in the A-term of the square lattice. The 0(7) invariance has 
some interesting implications for the spiral phases in a lightly doped antiferromagnet on the 
honeycomb lattice, as we will see in the next section. 

3. Applications 

In this section we will consider two nontrivial applications of the effective Lagrangian constructed 
above: The emergence of spiral phases in the staggered magnetization vector order parameter 
and the formation of two-hole bound states. 

3.1. Spiral phases 

Let us first investigate whether so-called spiral phases in the staggered magnetization can 
occur upon hole-doping. This problem has been studied before using microscopic theories 
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[33 S3 EH EH] . In Ref. [24] the problem was studied systematically in the case of a square 
lattice using effective Lagrangians. Here we consider the honeycomb lattice - details of the 
calculation can be found in Ref. }25j . 

In our study we assume that the holes are doped into the system homogeneously and thus 
limit ourselves to configurations of the staggered magnetization that are either homogeneous 
themselves or generate a constant background field Vi(x) for the charge carriers. As was 
shown in [24], the most general configuration of this kind represents a spiral in the staggered 
magnetization. We also assume that four-fermion contact interactions are weak and take them 
into account perturbatively - again, for details see Refs.[21j [24]. 

When one populates the various hole pockets, one must distinguish a total of four different 
cases. First we have a-pockets and /3-pockets. Moreover, the energy eigenstates of a hole pocket 



/ can acquire two values, E[_ and Ei_, where the former corresponds to the lower energy. One 
thus starts considering the case of populating all four hole pockets, i.e. with both flavors / = a, /3 
and with both energy indices ±. One then proceeds with populating only three, only two and, 
finally, only one hole pocket. This gives rise to various phases of the staggered magnetization, 
which can be either homogeneous or a spiral phase as depicted in figure 4. 
More precisely, the energy densities of the various phases take the form 



?f 



1 2 

€i = e + Mn + -Kin 



(24) 



Here eo is the energy density of the system at half-filling and n is the total density of holes. The 
index % refers to the number of hole pockets that are populated in the corresponding phase. The 
compressibilities «j are given by 
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Figure 5. The compressibilities Ki as functions of M' A 2 /2irp s determine the stability ranges of 
the various phases. A homogeneous phase, a spiral, or an inhomogeneous phase are energetically 
favorable, for large, intermediate, and small values of p s , respectively. 
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4irp s - M'A 2 1 r 

8(Gi + G 2 + G 3 )7rp s 



(3irp s - M'A 2 ) 2 16 
(4Gi + 3G 2 + 3G 3 )M'A 2 



«4 = ^ + \(G! + G 2 + G 3 ), (25) 

and shown in figure 5 as functions of M'A 2 /27rp s . For large values of p s , spiral phases cost a 
large amount of magnetic energy and the homogeneous phase is more stable. To be more precise, 
in this regime one has K4 < K3 < K2 < «i- Notice that k\ is always larger than K2 for any value 
of p s . As p s decreases and reaches the value 

M'A 2 (M') 2 A 2 G! (na . 
Ps = — + (26) 

at leading order in the 4-fermi couplings one finds K2 = = K4. For smaller values of p s , the 
two-pocket spiral is energetically favored until k 2 becomes negative and the system becomes 
unstable against the formation of spatial inhomogeneities of a yet undetermined type. 

Interestingly, unlike in the square lattice case, due to the accidental continuous 0(7) spatial 
rotation symmetry, at leading order a spiral does not have an a priori preferred spatial direction. 
It is instructive to compare the results presented here with the results obtained in the square 
lattice case |24j . Qualitatively the stability ranges of various phases are the same for both lattice 
geometries except that the one-pocket spiral is never energetically favored on the honeycomb 
lattice while it is favorable in a small parameter regime on the square lattice. 
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Figure 6. Tree-level Feynman diagram for one-magnon exchange between two holes. 



3.2. Two-hole bound states 

In the effective theory framework, at low energies, holes interact with each other via magnon 
exchange. Since the long-range dynamics is dominated by one-magnon exchange, we will 
calculate the one-magnon exchange potentials between two holes of the same flavor a and j3 
and of different flavor, and then address the question regarding the formation of two-hole bound 
states. The same problem was also considered in Refs. [401 Bdl |4"2] using microscopic theories. 

We expand in the magnon fluctuations m\{x) and m 2 {x) around the ordered staggered 
magnetization 

For the composite magnon fields this implies 

v^(x) = -^—d ll [m 2 (x)±im 1 (x)} + o(rn 3 ), 

Ps 

v^{x) = —— [mi(x)d^m 2 {x) - m 2 {x)d l j 1 mi(x)] +0 (m 4 ) . (28) 

Vertices with v^(x) involve at least two magnons, such that one-magnon exchange results from 
vertices with v^(x) only. As a consequence, two holes can exchange a single magnon only if 
they have anti-parallel spins (+ and — ), which are both flipped in the magnon-exchange process. 
We denote the momenta of the incoming and outgoing holes by p± and p±, respectively. The 
momentum carried by the exchanged magnon is denoted by q. The incoming and outgoing holes 
are asymptotic free particles with momentum p = (pi,p 2 ) and energy E(p) = M + pf /2M' . 
One-magnon exchange between two holes is associated with the Feynman diagram in figure 
6. Evaluating these Feynman diagrams, one finds the resulting potentials for the various 
combinations of flavors /, / G {a, /?} and couplings F, F G {A, K}. 

The leading contribution to the low-energy physics, as we noted earlier, comes from the A 
vertex. From here on, we therefore concentrate on the potential with two A vertices only. In 
coordinate space the AA-potentials are given by 

(rirl\vH\r + r„) = v£(r) 5(r + - fl) «5(r_ - r±), (29) 

with 



V[[{r) = -§-/ 2) (n, V[[{r) = _^exp(2 J a / ^). (30) 



Here r = f + — r_ denotes the distance vector between the two holes and (p is the angle between 
r and the xi-axis. The <5-functions in Eq. (|29p ensure that the holes do not change their position 



during the magnon exchange. Note that the one-magnon exchange potentials are instantaneous 
although magnons travel with the finite speed c. Retardation effects occur only at higher orders. 

Interestingly, in the AA channel, one-magnon exchange over long distances between two 
holes can only happen for holes of opposite flavor. For two holes of the same flavor, one-magnon 
exchange acts as a contact interaction. Here we concentrate on the long-range physics of weakly 
bound states of holes and therefore will only consider the binding of holes of different flavor. 

Let us therefore investigate the Schrodinger equation for the relative motion of two holes with 
flavors a and (3. In the following, we will take care of short distance interactions by imposing a 
hard-core boundary condition on the pair's wave function. Due to the accidental Galilean boost 
invariance, without loss of generality, we can consider the hole pair in its rest frame. The total 
kinetic energy of the two holes is given by 

f=a,P f=a,P 

We now introduce the two probability amplitudes ^i(r) and ^(j) which represent the two 
flavor-spin combinations a+/3_ and a_/3+, respectively, with the distance vector r to point from 
the f3 to the a hole. Because the holes undergo a spin flip during the magnon exchange, the two 
probability amplitudes are coupled through the magnon exchange potentials and the Schrodinger 
equation describing the relative motion of the hole pair is a two-component equation. Using the 
explicit form of the potentials, it takes the form 

( -t^A 7 ^exp(-2^) \ / vD^f) \ _ / ^(f) \ 
^ 7^ expQtop) A ) { Mr) )~ { * 2 (r ) ) ' 

with 

A 2 

1=~ • (33) 

2irp s 

As it turns out, magnon-mediated forces can lead to bound states only if the low-energy constant 
A is larger than a critical value given by 



A, - ^ (-) 

Interestingly, the same critical value also arose in the investigation of spiral phases in a lightly 
doped antiferromagnet on the honeycomb lattice in Eq. (|26p . There it marked the point where 
spiral phases become energetically favorable compared to the homogeneous phase. 
With the separation ansatz 

*l(r, <p) = Ri(r)exp(imi(p), *2(f, <p) = R%{r) exp(im2f), (35) 

the relevant equation which is associated with an attractive potential and thus potentially leads 
to the formation of bound states amounts to 



d 2 1 d \ M lWN 1 



R(r) = -M'\E\R(r), (36) 



with E = —\E\ and R{r) = R\{r) — R<z{r). It has to be pointed out that this equation refers to 
the case m\ = — 1 and mi = 1 where the system (|3"2"j) can be decoupled (see [2"2"]). 

The same equation also occurred in the square lattice case [20l [23] and can be solved along 
the same lines. As it stands, the equation is ill-defined because the 1/r 2 potential is too singular 



at the origin. However, we have not yet included the contact interaction proportional to the 
4-fermion coupling G3. Here, in order to keep the calculation analytically feasible, we model the 
short-range repulsion by a hard core radius Tq, i.e. we require R(tq) = for r < rn. Eq. (|36p is 
solved by a modified Bessel function 



R(r) = AK v (yJM'\E\r), v = i^yM' - 1, (37) 

with A being a normalization constant. Demanding that the wave function vanishes at the hard 
core radius gives a quantization condition for the bound state energy. The quantum number n 
then labels the n-th excited state. For large n, the binding energy is given by 

1 / -2vrn \ . , 

E n ~ -T7T2 ex P / tt, \ • ( 38 ) 



M'rl VV7M 7 - 1 

The binding is exponentially small in n and there are infinitely many bound states. While 
the highly excited states have exponentially small energy, for sufficiently small vq the ground 
state could have a small size and be strongly bound. However, for short-distance physics the 
effective theory should not be trusted quantitatively. Still, as long as the binding energy is 
small compared to the relevant high-energy scales, our result is valid and receives only small 
corrections from higher-order effects such as two-magnon exchange. 

We now turn to the discussion of the angular part of the wave equation. The ansatz (|35|) 
leads to the following solution for the ground state wave function 

•<™> = ( %% = ) • (») 

Under the 60 degrees rotation O, using the transformation rules of Eq. (|17|) . one obtains 

°y(r,<p) = -V(r,<p). (40) 

Interestingly, the wave function for the ground state of two holes of flavors a and j3 thus exhibits 
/-wave symmetry^] The corresponding probability distribution depicted in figure 7, on the 
other hand, seems to show s-wave symmetry. However, the relevant phase information is not 
visible in this picture, because only the probability density is shown. Interestingly, for two-hole 
bound states on the square lattice, the wave function for the the ground state of two holes 
of flavors a and f3 shows p-wave symmetry, while the corresponding probability distribution 
(which again does not contain the relevant phase information) resembles d x i_ y 2 symmetry |20j . 
Remarkably, the ground state wave function (139[) of a bound hole pair on the honeycomb lattice 
remains invariant under the reflection symmetry R, the shift symmetries Di, as well as under 
the accidental continuous rotation symmetry 0(7). 

We find it quite remarkable that the /-wave character of the two-hole bound state on the 
honeycomb lattice is an immediate consequence of our systematic effective field theory analysis. 
The question regarding the true symmetry of the pairing state realized in the dehydrated version 
of Na2Co02 x y H2O, still seems to be controversial [33]. Still, a careful analysis of the available 
experimental data for this compound suggests that the pairing symmetry indeed is /-wave [33]. 

4. Conclusions 

The effective theory for the insulating antiferromagnetic precursors of high-temperature 
superconductors is the condensed matter analog of baryon chiral perturbation theory: Magnons 



1 Strictly speaking, the continuum classification scheme of angular momentum eigenstates does not apply here, 
since we are not dealing with a continuous rotation symmetry. 



Figure 7. Probability distribution for the ground state of two holes of flavors a and f3. 



play the role of the pseudoscalar mesons, while the holes and electrons are the analog of the 
baryon octet. We have analyzed the symmetries of the underlying Hubbard and t— J-type models 
in detail and reviewed the systematic construction of the effective low-energy field theory for 
weakly hole-doped antiferromagnets on the honeycomb lattice. The procedure is fully systematic 
order by order in a derivative expansion of the magnon and hole fields. As two nontrivial 
applications we have considered the existence of spiral phases in the staggered magnetization 
order parameter and the formation of two-hole bound states mediated by magnon exchange. 

Regarding spiral phases, we have limited ourselves to constant composite vector fields Vi(x)' 
which implies that the fermions experience a constant background field. Unlike in the square 
lattice case, due to the accidental continuous 0(7) spatial rotation symmetry, at leading order a 
spiral does not have an a priori preferred spatial direction. However, since the 0(7) symmetry 
is broken explicitly by the higher-order terms, once such terms are included, one expects the 
spiral to align with a lattice direction. We also investigated the stability of spiral phases in the 
presence of 4-fermion couplings. If these couplings can be treated perturbatively, for sufficiently 
large values of p s , the homogeneous phase is energetically favored. With decreasing p s , a two- 
pocket spiral becomes energetically more favorable. In contrast to the square lattice case, the 
one-pocket spiral is never favored. For small values of p s the two-pocket spiral becomes unstable 
against the formation of inhomogeneities of a yet undetermined type. 

Regarding the formation of two-hole bound states mediated by magnon-exchange, we have 
studied the effect of the magnon-hole vertex. Again in contrast to the square lattice case, it 
turned out that the magnon-hole coupling constant A must exceed a critical value in order to 
obtain two-hole bound states. Our analysis implies that the wave function for the ground state 
of two holes of flavors a and (3 exhibits /-wave symmetry (while the corresponding probability 
distribution seems to suggest s-wave symmetry). This is quite different from the square lattice 
case, where the wave function for the ground state of two holes of flavors a and (3 exhibits p-wave 
symmetry (while the corresponding probability distribution resembles d x 2_ y 2 symmetry). 

It is quite remarkable that all these results unambiguously follow from the very few basic 
assumptions of the systematic low-energy effective theory, such as symmetry, locality, and 
unitarity. The effective theory provides a theoretical framework in which the low-energy 
dynamics of lightly hole-doped antiferromagnets can be investigated in a systematic manner. 
In particular, after the low-energy parameters have been adjusted appropriately, the resulting 
low-energy physics is completely equivalent to the one of the Hubbard or t-J model. 
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